The making of ferromagnetic Fe doped ZnO nano-clusters 
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In this letter, the authors present a study of the energetics and magnetic interactions in Fe doped 
ZnO clusters by ab-imtio density functional calculations. The results indicate that defects under 
suitable conditions can induce ferromagnetic interactions between the dopant Fe atoms whereas 
antiferromagnetic coupling dominates in a neutral defect-free cluster. The calculations also reveal 
an unusual ionic state of the dopant Fe atom residing at the surface of the cluster, a feature that is 
important to render the cluster ferromagnetic. 
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Mn doped GaAs is a prototype diluted magnetic semi- 
conductor (DMS). However its Curie temperature could 
be raised only to a maximum of ^^170 making it 
hardly useful for room temperature spintronic appli- 
cations. For room temperature ferromagnetism (FM) 
ZnO and GaN were predicted to be promising as host 
materials.^ Following this theoretical prediction there has 
been an avalanche of experimental work to realize room 
temperature FM in these systems. Among these mate- 
rials transition metal (TM) doped ZnO received special 
attention as it was possible to dope TM's not only in bulk 
ZnO but also in thin films as well as nano-crystalline form 
of ZnO. However the experimental results on the realiza- 
tion of room temperature FM are highly controversial 
but there is general consensus that FM strongly depends 
on methods and conditions employed in the preparation 
of the samples. 

Nanoclusters of ZnO capable of sustaining ferromag- 
netic order are of particular interest since a clear under- 
standing of the finite size effect on the magnetic proper- 
ties of such systems is essential for the development of 
high-density storage media with nano sized constituent 
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FIG. 1: (Color online) (a); Formation energy of one Fe atom 
doped at a Zn site of a ZnO cluster as a function of excess 
chemical potential of O {^i'^), (b): Energy difference (Ai5) 
between antiparallel (AFM) and parallel (FM) magnetic ori- 
entations when two Fe atoms are doped at the Zn sites as 
a function of various configurations (with increasing distance 
between two Fe ions in the cluster). Inset: Relaxed structure 
of the pure Zni20i2 cluster. 



particles or crystallites. In contrast to bulk the nano par- 
ticles are dominated by surfaces and are therefore nat- 
urally susceptible to defects that may have crucial im- 
pact on their magnetic properties and future application 
in a spintronic device. Kittilstved et al.^ have shown 
the importance of defects for ferromagnetism in Mn and 
Co-doped ZnO nano-crystals. Very recently, Karmakar 
et al* have indicated that the coexistence of Fe^+ and 
Fe'^"'' states promotes ferromagnetism in Fe doped ZnO 
nano-cluster. They argued that the surface electronic 
structure of these nanoparticles is quite different from 
the core presumably due to defects and is dominated by 
Fe'^"'" ions, which is crucial for the magnetic properties. 
These observations are further supported by recent x-ray 
magnetic circular dichroism (XMCD) studies on the same 
samples.^ In the present letter, using ab-initio electronic 
structure calculations we have investigated the conditions 
that will render Fe doped ZnO clusters to be ferromag- 
netic. 

For the electronic structure calculations we have used 
the projector augmented wave (PAW) methodS as im- 
plemented in the VASP package.^' - PAW potentials with 
12 valence electrons for Zn {Sd^^As"^), 6 for O {2p^2s'^) 
and 8 for Fe {3d^As^) with an energy cutoff of 500 eV for 
the plane wave expansion of the PAW's was employed in 
our calculations. The exchange-correlation part was ap- 
proximated in the generalized gradient approximation^ 
including a Hubbard U for the dopant Fe. For Fe, U (on- 
site d — d Coulomb interaction) — 4.0 eV and J (on-site 
exchange interaction) = 1.0 eV was chosen. The forma- 
tion energy (FE) of the defects was calculated using the 
following expressionjiSiiiii^ 

FE - S(Zn„,0„(a,(7)) -S(ZnpOp(pure)) 

-I- ^Ualla + qiEy + ep) (1) 
a 

where a is the defect atom added or removed from 
the pristine ZnO, ria is the number of each defect 
atoms: Ua — —1 (+1) for adding (removing) one atom, 
i?(ZnniOn(a, g)) is the total energy with defect a and 
charge q, while £'(ZnpOp(pure)) is the total energy of the 
pure ZnO cluster, ep is the Fermi level measured with re- 
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spect to the energy of the highest occupied molecu- 
lar orbital (HOMO) in the pristine ZnO cluster. /Xq is the 
chemical potential of atom a in some suitable reservoir. 
Chemical potentials of Zn and O are obtained from the 
condition: p(^|„ + + Mo + Mo) = -B(ZnpOp(purc)), 
where and iiq are the energies per atom of Zn and 
O respectively in the standard elemental structure and 
/if „ and Hq are excess chemical potentials of Zn and O 
which obey the above constraint at equilibrium. In ad- 
dition, the excess chemical potentials should satisfy the 
relation /xf „ < and < to assure that Zn and O 
in the elemental structure are not formed. Depending on 
growth condition, the excess chemical potentials can have 
two extreme limits either /if „ — oi — referred to 
as the Zn rich and O rich limits respectively. For Fe 
we have taken the chemical potential to be the energy 
per atom in the most commonly occurring bcc structure 
and assumed /i|,g — 0. The total energies of the charged 
simulation cell were computed by compensating any ad- 
ditional charge by a uniform jellium background. 

In view of the recent reportai^^ii^ that Zni20i2 is 
the magic cluster with the most stable configuration we 
confined our study to the 24-atom Zni20i2 cluster. The 
cluster was simulated in a large cell (a cube of 20 A) sur- 
rounded by vacuum and periodic boundary conditions. 
In order to obtain the optimized structure we considered 
different initial structures and performed structural re- 
laxations which led to the optimized structure displayed 
in Fig. l(see Inset). The structure is spherical with all the 
atoms on the surface of the sphere providing an unique 
opportunity to study the crucial surface effects. All the 
atoms are very symmetrically arranged with six 4 atom 
and eight 6 atom rings. The diameter of the cluster is 
~6.35 A. The HOMO-LUMO gap for this cluster was 
found to be 2.34 eV. 

One atom of the cluster was substituted with an Fe 
atom and this cluster was fully relaxed. We find that 
substitution of Zn atom by Fe is energetically more fa- 
vorable than substitution of an O atom in the ZnO clus- 
ter. The formation energy of doping one Fe atom at the 
Zn site in the cluster has been calculated and is plotted 
in Fig. 1(a) as a function of excess chemical potential of 
Oxygen (/i^). As expected, the formation energy of dop- 
ing Fe at the Zn site is higher in the Zn-rich limit. The 
singly doped cluster is magnetic with a magnetic moment 
of 4 fiB, as expected from Hund's rule. 

In order to check the tendency of magnetism (ferro- 
magnetic (FM) or antiferromagnetic (AFM)) a pair of 
Fe ions are substituted into the cluster assuming a paral- 
lel (FM) as well as an anti-parallel (AFM) arrangement 
of magnetic moments and the structure of these clusters 
are relaxed. The difference in total energies between the 
AFM and FM arrangement of magnetic moments (Ai?) 
is displayed in Fig. 1(b) for seven possible configurations 
with increasing separation between Fe ions in the cluster. 
This energy is a measure of the inter-atomic exchange in- 
teraction. We observe that the coupling is always AFM 
and decays very fast with distance. This short-ranged 
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FIG. 2; (Color online) (a): Formation energies for one Fe 
atom doped at a Zn site in the cluster in the oxygen rich 
limit with and without vacancies for different charged states. 
Only the lowest formation energy of a particular charged state 
is shown. When a vacancy is present, the Fe ion and the va- 
cancy are at the closest possible positions. The dotted verti- 
cal line indicates the calculated band gap for the cluster; (b), 
(c), (d): corresponds to AE between AFM and FM magnetic 
orientation when a pair of Fe atoms at the closest possible 
position are doped at the Zn sites in presence of no vacancy, 
Zn vacancy, and O vacancy respectively. 
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FIG. 3: (Color online) (a): Partial DOS of Fe-d states and 
(b): charge density in an energy window corresponding to one 
minority spin state (indicated by a dashed arrow) in a cluster 
where a pair of Fe atoms are doped at Zn sites near a Zn 
vacancy with charged state —1. 



AFM superexchange interaction has also been observed 
in the case of transition metal doped bulk ZnOi^ This is 
due to the presence of transition metal states in the en- 
ergy gap of the host cluster. Our calculations reveal that 
it is not possible to stabilize ferromagnetism by doping 
Fe in a neutral charge state in pristine ZnO. However the 
abundance of room-temperature ferromagnetism in Mn, 
Fe and Co doped ZnO nanocrystals suggests that defects 
may play a crucial role to stabilize ferromagnetism in 
these systems. 

Therefore the formation energies were calculated for a 
single Fe atom doped into the cluster with and without 
defects (Zn and O vacancy) assumed in different charged 
states. We considered both Zn rich and O rich condi- 
tions. The Fermi energy is varied from the HOMO level 
of the pristine ZnO cluster to the value of the gap of bulk 
ZnO. In Fig. 2(a) we display our results for the formation 
energy of Fe doped ZnO under oxygen rich conditions. 
From this figure we note that near the p-type region. 
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Fe doped ZnO (without defects but in the charged state 
+1) has the lowest formation energy while in the n-type 
region, the lowest formation energy corresponds to the 
doped system with Zn vacancy in the charged state — 1. 
The stability of the above two configurations provide an 
important evidence that the ionic state of the dopant Fe 
atom residing on the surface of the cluster may be Fe'^+ in 
agreement with recent experimental results^ and may be 
attributed to the different co-ordination structure of the 
surface Fe atoms in comparison to the bulk terminated 
ones due to loss of ligands and co-ordination defects aris- 
ing from lattice distortion. In the Zn-rich limit the in- 
clusion of defects does not lower the formation energy in 
either n-type or p-type region. 

We have now examined whether the above two sta- 
ble configurations may support the tendency for ferro- 
magnetism. In Fig. 2(b),(c),(d) we have displayed the 
energy difference between the AFM and FM configura- 
tion by including two Fe atoms at the closest possible 
positions without and with Zn and O vacancies in vari- 
ous charged states. We gather from the figure that for 
both of the above mentioned stable configurations there 
is a tendency for ferromagnetism (See Fig. 2(b), (c)). The 
magnetic moment resides in each case predominantly on 
the dopant Fe atoms spilling over a bit onto the neigh- 
boring oxygen atoms. A pair of Fe atoms substituting Zn 
in the charged state —1 and -f 1 in the presence and ab- 
sence of a Zn vacancy respectively point to the fact that 
the dopant Fe atoms are in the ionic states Fe^+ {(fi) and 
Fe'^"'' (d^) (mixed valent). This ionic configuration is im- 
portant to stabilize ferromagnetism as was speculated in 
a recent experiment This is further illustrated in Fig. 3 
where we have considered a particular case namely a pair 
of Fe atoms substituted into the ZnO cluster in the pres- 
ence of a Zn vacancy in the charged state —1. For this 



case we plotted the Fe-d partial density of states (DOS) 
and show the charge density in a small energy window 
corresponding to a single Fe-d state in the minority spin 
channel as indicated in the figure. From the DOS we see 
that there is an admixture of the Fe-d states with the 
0-p states. Furthermore in the majority spin channel 
both the Fe atoms are completely occupied whereas in 
the minority spin channel only one Fe d state is occupied 
resulting in a net magnetic moment of 9 hb- The corre- 
sponding charge density for this Fe-d state in the minor- 
ity spin channel reflects a hopping induced hybridization 
between the neighboring Fe ions which is responsible for 
promoting ferromagnetism by a reduction in energy in 
the process of hopping. Although the oxygen vacancies 
do not have the lowest formation energies, their effect 
on the magnetic properties is also very interesting as it 
induces ferromagnetic interactions in otherwise antifer- 
romagnetically coupled magnetic atoms in a defect-free 
system (See Fig. 2(d)). 

In conclusion, we have studied the energetics and mag- 
netic interactions in an Fe doped ZnO nano-cluster. It 
has been found that in a pure ZnO cluster the Fe atoms 
couple antiferromagnetically. The presence of Zn and 
O vacancies induces intra-cluster ferromagnetic coupling 
and hence a net magnetic moment for the cluster. We 
have argued that the ionic state of the dopant Fe ion 
plays a crucial role to make the cluster ferromagnetic. 
Thus one can envisage the engineering of defects in nan- 
oclusters to design spintronic devices. 
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